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Abstract 

We constrain the thermal evolution of the universe with a decaying cosmological term by using 
the method of the analysis for the Wilkinson Microwave Anisotropy Probe (WMAP) observation 
data. The cosmological term is assumed to be a function of the scale factor that increases toward 
the early universe, and the radiation energy density is lower compared to that in the model with the 
standard cosmological constant (ACDM). The decrease in the radiation density affects the thermal 
history of the universe; e.g. the photon decoupling occurs at higher-z compared to the case of 
the standard ACDM model. As a consequence, a decaying cosmological term affects the cosmic 
microwave background (CMB) anisotropy. Thanks to the Markov-Chain Monte Carlo method, we 
compare the angular power spectrum in the decaying ACDM model with the CMB data, and we 
get severe constraints on parameters of the model. 

PACS numbers: 95.36.+X, 98.80.-k, 98.80.Es 
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I. INTRODUCTION 

Recent astronomical observations such as high redshift type la supernovae (SNIa) 
cosmic microwave background (CMB) anisotropy [s^, suggest the existence of dark energy 
strongly. Although many researchers have investigated dark energy, its nature is still un- 
known. In the proposed n^ode. „. the equation ot state ot da* energy g the cosntological 
constant cannot be excluded [2|]. 

If we assume that the dark energy is equivalent to a cosmological constant A, there rises 
again so called a cosmological constant problem^: the present value of A is extraordinarily 
small compared with an inferred vacuum energy during the Planck time. To solve this 
problem, it is natural to consider that A decreases from a large value at the early epoch to 
the present value. Many functional forms of A have been suggested : for instance decaying-A 
las been introduced as a function of a scalar field in the Brans-Dicke gravitational theory 
6j. The evolution of the universe under various models of A which are included in the 
energy-momentum tensor has been investigated analytically |7|. 

In addition, interacting A with other kinds of energy has been also discussed. The vacuum 
energy of A coupled with baryon could be ruled out, because baryon-antibaryon created 
through vacuum decay causes pair-annihilation. The produced high energy gamma ray flux 
is contradicted with the observations [8]. On the other hand, the vacuum energy decayed 
into the photon could affect the cosmological evolution significantly. Assuming that the ratio 
of the vacuum energy to the radiation is constant at the radiation dominated era {z > 10^), 
Freese et al. js] investigated the effects on primordial nucleosynthesis and obtained a limit of 
vacuum to photon energy ratio, which is less than 0.07. Furthermore, it is pointed out that 
observational constraints from the CMB intensity put the limit on the ratio of the vacuum 
to the radiation energy to be ~ 10~^ 

From the thermodynamical point of view, the temperature-redshift relation is modified 



by including adiabatic photon creation due to vacuum decay [10 



A phenomenological 



decaying-A has been found to affect the cosmological evolution after the recombination 



12 



13 



14( 1 . In models having A terms as a function of the scale factor, the radiation and 



matter temperatures would be significantly lower compared to the standard cold dark matter 
model with a constant A (SACDM) [12]: the molecular formation is occurred at earlier epoch 
by Az < 10^ [3]. Furthermore, it is shown that in some parameter regions, the radiation 
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temperature could become higher compared with the SACDM model, which is found to be 

n 

consistent with the observational result of z < 4 |14j |. 

Related to the recent observations, the first star formation that occurred at the end 

n 

of the dark age has been investigated progressively [15]. In decaying A models, using a 
cooling diagram, the first star formation was estimated to occur at an earlier epoch by 
Az ~ 20 with its mass ~ IO^Mq |13|]. In the meantime, from an observational approach, 
the CMB polarization observed by the WMAP satellite predicts via measured reionization 
redshift with use of the SACDM model that a first object was formed around z = 10 
[sl, 3]. Since the CMB anisotropies give severe constraints on parameter regions concerning 
the cosmological evolution, we can also estimate the era of the first star formation in the 
decaying-A cosmology using the CMB anisotropy. 

In the present paper, we constrain the parameter regions that determine the thermal 
history of the universe with a A decaying into the photon (hereafter we call it DACDM 
). In Sec. nil we describe briefly the thermal evolution and clarify the effects on photon 
decoupling in the DACDM model. In Sec. Illll we examine the consistency between DACDM 
and the CMB anisotropy data of Winlinson Microwave Anisotropy Probe (WMAP) using 
the Markov-Chain Monte Carlo (MCMC) method. Summary and discussion are given in 
Sec. Ml 



II. THERMAL EVOLUTION WITH A DECAYING COSMOLOGICAL TERM 

Using the Friedmann- Robertson- Walker metric, the Einstein equation and/or the energy- 
momentum conservation law are written as follows: 

[a) = W 

f, = -3^{p+p), (2) 

where a, K and G are the cosmic scale factor, the curvature and the gravitational constant, 
respectively. We note that bars such as p and p indicate the average values during the 
cosmological evolution. We choose the unit such that the velocity of light c = 1. Note that 
dots in Eqs. ([1]) and ([2]) indicate the derivative concerning a conformal time r. The total 
energy density p and the pressure p are written as 

P = Pm + P^+Pu + PA, P = P^+Pu+ Pa, (3) 
3 



where the subscripts m, 7, i/, and A indicate the nonrelativistic matter (baryon plus cold 
dark matter), photon, neutrino, and a cosmological term, respectively. The equation of 
states p/ p for individual components are written as. 



v/p 



1/3 relativistic particles, 

non-relativistic particles, 

— 1 cosmological term . 



(4) 



Here the energy densities of matter and neutrino vary as = Amo*^"^ ^"^^ Pu — Puo^~'^j 
where the subscript means the present value. 

From Eqs. ([2]), ([3]) and (HI), we get the evolution equation of the photon energy density 
after the epoch of electron-positron pair-annihilation : 

a da 

with the density parameter fi- 

Pi 



dn n 

da 



(5) 



p 



3Hi 



crit 



Pcrit ■ SttG 

where i^Q is the Hubble constant in units of km/sec/Mpc. 

In the DACDM model, the evolution of the photon is affected by the time-dependent 



cosmological term. In this work, we assume a functional form of A as follows [l2, 13, Q, 'v\ : 



(6) 



where ^^^^,^^2 m are constants. Note that the present value of fi^ is expressed by 



^AO ~ ^ 



-r ^2^2- Formalism of this paper is based on previous studies in Refs. 
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Integration of Eq. ([5]) with IQ and Stefan-Boltzmann's law, oc T^, leads to the following 
photon temperature as a function of the scale factor jl^ . 
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where is the present photon temperature and a = ■mQ/<^2/{4: — rn). The present photon 
energy density Q^q = 2.471 x 10^^ h^'^iT^^/ 2. 725 K)^ with the normalized Hubble constant h 
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FIG. 1: Upper panel: the evolution of the photon temperature in DACDM (m = — 1.2) with 
— 10~^ after the hydrogen recombination era. Lower panel: the ratios of T-y to that of SACDM. 

{Hq = 100 h km/sec/Mpc). Second terms in Eq. (JTj) are characteristic ones in the DACDM 
model; the evolution of the photon is modified by those terms. 

Figure [1] illustrates the evolution of the photon temperature in the DACDM model. It 
can be seen that Ty in DACDM is lower compared to that in SACDM. For < m < 4, the 
photon evolves as oc at the early epoch and the slope of against a decreases due 
to the contribution of a^"™" near the present epoch. For m > 4, the opposite results are 
obtained. Therefore, decaying-A affects the thermal history at around the present epoch of 
z < 10^ and particularly after the hydrogen recombination era. 

If and/or m is very large, the solution with present T^q fixed indicates that the 
total energy density becomes negative for some epoch of z > 0. As a consequence, these 
parameters were constrained as mfiA2 < 10"'^ [ij]. In our analysis, the photon temperature 
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for large m and/or becomes also negative at some epoch of a < 1. By excluding this 
kind of solution, we obtain the upper limits on both $7^2 ^'^d m from Eq. ([7]): 

a < (rn < 4). (8) 

In the case of m > 4, we assume > until the primordial nucleosynthesis epoch, a = 
10~^°, and obtain the limits 

fi,o > 920^2, (rn = 4) (9) 

fi^o > -10'°('"-^)a (m > 4) (10) 

On the other hand, for i7^2 < or m < 0, we find that becomes negative at some time 
of a > 1. Therefore we impose the two conditions of ^2^2 — ^"^^ m > 0. 

In the DACDM model, the cosmological term decreases from the early time to the present, 
because the second term in Eq. ([6]) becomes larger than the first one at the early epoch of 
z > 0. Since the first term in Eq. ([6]) dominates near the present epoch, the A term is 
nearly constant for low- 2;. Although cosmological models with the A term have been tightly 
constrained from the luminosity-redshift relation of SNIa, effects on the expansion rate are 
neghgible in the DACDM model. 

In the SACDM model, the ratio of the photon to the neutrino temperature is T^/T^ = 
(11/4)^/^ after electron-positron annihilation, because and evolve as oc a^^. In the 
DACDM model, since a decaying-A alters the evolution of the photon, the ratio of the 
photon-to-neutrino depends on time. Nonetheless if we set T^/T^, = (11/4)^/^ at a = 10^^° 
in the DACDM model, the present neutrino temperature is lower than that in SACDM as 
seen in the second column of Table [H Although recent observational studies put constraints 
on properties of the cosmic neutrino background such as the neutrino species or masses (e.g. 
jsl, 13])' there is no observation about its temperature (or the energy density). Therefore, 
DACDM seems to have no problems for the lower neutrino temperature. 

CMB temperature at z = was measured accurately by the Far Infrared Absolute Spec- 
trophotometer of Cosmic Background Explorer: T = 2.725 ± 0.002 K at 2cr C.L. [l9|. On 
the other hand, CMB temperature observation at 2; > are reported in literature [2C 



Consistency of the temperature evolution in the DACDM model with these observational 
results has been discussed [l^. When m and/or 17^2 take a large value, in the DACDM 
model is not consistent with the temperature observation. Puy IJ] has put constraints only 



on the m — Q 



plane from the temperature observation of 2; < 1 2l| : 



m| < 1, |fi 



< 10 



-4 



(11) 



'A2I ^ 



These hmits are obtained by comparing the observational temperature included la error 
with the temperature given in Eq. ([7j). Thus, we perform the extended analysis to constrain 
the parameters using the available observations as precisely as possible. Figure [2] shows 
constraints on the m — Q^2 plane from the observational temperature using the same analysis 
in Ref. under the theoretical request of > 0. Constraints from Tj at z > 1 are 
similar for m > 4 as shown in the upper panel of Figj2] and those from at z < 1 have 
large uncertainty as shown in the lower panel of Figj2l These limits of m or i7^2 obtained 
from the observed temperatures are consistent with the excluded region by Eqs.®, (j9]) and 
f ITU]) . Since obtained constraints are rather rough, we put further severe limits using CMB 
anisotropy as shown in Sec JIIIi 

Now we estimate the epoch of photon decoupling by calculating the visibility function 
g{t) which has a peak at the epoch of the last-scattering surface: 



where ctt is the Thomson-scattering cross-section and is the number density of the free 
electrons that depends on the recombination history of the universe. 

Figure [3] illustrates the visibility function as a function of redshift in the SACDM and 
DACDM models. The epoch of photon decoupling shifts to higher-z as shown in the third 
column of Table [B When we take the upper limits obtained from Eq. (IH]), (fiA2,m) = 
(10"'*, 1.2), the photon decoupling occurs at Zdec = 2040 that is earlier by Azdec ~ 950 
compared to the case of the SACDM. In addition, we find that photon temperature at the 
last-scattering surface is about 0.1% lower than that in the SACDM model, as shown in the 
fourth column of Table [H These effects should be further constrained by the cosmological 
observations such as CMB anisotropy as shown in the following sections. 

III. EFFECTS ON CMB ANISOTROPY AND CONSTRAINTS BY MARKOV- 
CHAIN MONTE CARLO ANALYSIS 

Before discussing effects of a decaying-A on the CMB power spectrum, let us formulate 
the Boltzmann equation for the photon in the DACDM model based on the cosmological 
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FIG. 2: Constraints on the m — plane from observational temperatures. The black-solid line 
shows the upper limits of parameters by Eqs. ([8]), Q, and (jlOp . and the other lines indicate 
the upper limits obtained from the same analysis in Ref. jl4| . Upper panel: constraints from the 
temperature at z > 1. Lower panel: constraints from the temperature at z < 1. 
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FIG. 3: Visibility function vs the redshift for a fixed parameter ^^^2 = 10~^. Peaks of g show the 
epoch of photon decouphng. As the values of the parameter m increase, the photon decoupling 
occurs at an earlier epoch. 

perturbation theory. The line element in the synchronous gauge with the flat space is written 

as 

ds"^ = a{T) [-dr'^ + {6ij + hij) dx'dx^] , (12) 

where h^j is the metric perturbation. We introduce two fields h{k, r) and ri{k, r) in the 
Fourier /c-space and write the scalar mode of h^^ as the Fourier integral: 

5. 



hij{x, t) = d ke 
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TABLE I: Present neutrino temperature {z = 0), redshift, and photon temperature at the last- 
scattermg surface {z^ec) in the decaying-A model with Oa2 = 10"'*. 



Parameter 



Tu [K] {z = 0) 



^dec 



[K] (Zde 



m = 0.0 
m = 0.5 
m = 1.0 
m = 1.2 



1.945 
1.774 
1.416 
1.022 



1087 
1188 
1480 
2043 



2965 
2957 
2939 
2921 



where kj = kk- with the unit vector k- . Components of the energy-momentum tensor with 
perturbed parts are given by 

T° = -(p + 5p), 

T^ = -n = {p + P)v,, (13) 
T^ = {p + 5p)5] + T}., 

where 5p, Sp, Vi and are perturbed parts of the energy density, pressure, velocity of fluids, 
and anisotropic stress, respectively. 

In the DACDM model, the energy density of the photon in the background part is ob- 
tained from Eq. ([5]). The perturbed parts of the equation of the energy-momentum conser- 
vation that correspond to the first order perturbation reduce to the following equations, 

9, = \k\ - ea-, + ^9^, (15) 
where = Sp^/p^ and 6'^ is the divergence of the fluid velocity, 6*^ = ik^vj. is defined by 

[p-y + Pj) 0-7 = - (^k ■ kj - ^6ij^ S}. 

Equations ([H]) and (fT5l) are the continuity and Euler equations, respectively. Note that we 
take into account the interaction terms proportional to pa in these equations and the effects 
will be small since the ratio P/^/p-y at z = 10^ is less than 10~^. 

Moreover, to construct the perturbed evolution equation of the photon, we need the 
contribution of the higher multipole moments. The Boltzmann equation for a relativistic 
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particle in /c-space is written as follows 22| : 



, ,a/o dlnL { . h + 6T] A fdf\ 



col 



where ^ = k ■ fi, and = gn, is the 3-dimensional momentum. The right-hand side of 
Eq. fll6l) is the collision term. The distribution function expressed by the convolution of 
zeroth-order and the perturbed part is written as 

f{x\q,nj,T) = /o(g,r)(l + *(x\g,nj-,r)). 

In the SACDM model, the zeroth-order distribution function, /q, of the photon is described 
as f^Q^q) = l/(/ipexp (g/fc^T^) — 1) , where and are the Planck and Boltzmann 
constants, respectively. However, in models with a created photon by decaying vacuum 
energy, the spectral distribution of CMB is a function of both temperature and the comoving 
number of photon, where /o takes a generalized Planckian form [11]. Therefore, we cannot 
drop the second term in the left-hand side in Eq. (fT6|) . which is the time derivative of /q. 

To obtain the Boltzmann equation for photon, we expand the angular dependent part of 
the perturbation in a series of Legendre polynomials P;(fc ■ n) as follows: 

F,{k, n, r) ^ Wo(g,^)^ = f2 i-^y (2^ + 1) F.ii^, r)PM. 
^ iq^dqf^[q,T) ^ ^ 

We integrate Eq. fll6p multiplied by q^dqf^ over the whole p space and divide it by J q^dqf'Q. 
Then we obtain the following Boltzmann equation for the CMB photon in A;-space: 

+ ikf,F^ + ^ + 6^) ^2(P) + lh- ^F^ = (^^^ ^ . (17) 

The last term of the left-hand side in Eq. f|T7|) corresponds to a new one that appeared 
in the DACDM model. The collision term in the right-hand side is described as Thomson 
scattering {2^ : 



dr 



+ F^^ + Ah-v^-]^ {F^ + G^o + ^72) A(P) 



:i8) 



col 

where G^i is the difference of the two linear polarization components. 

Substituting the Legendre expansion for F^, and using the orthonormality of the Legendre 
polynomial with the recursion relation (/ + l)P^_^]^(yu) = {21 + l)fiPi{fi) — lPj^_^{n) , we get 
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the Boltzmann equations for the photon in the DACDM model as follows: 



1 



t P7 



(19) 



a. 



7 



15 ^ 

F = 
where 
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Pa 



10 



20 



~ 10^^'^^ + YS + ^ 77ia^ea;g^yf^7 + ^^Ti^'^e^e^T (^70 + '^72) + —^'^^) 



(/ + 1) - an.x.a^F^, + ^F^, (/ > 3), 

P7 



(21) 



-^70; ^7 



4^^7i' S 



We calculate the CMB power spectrum by modifying the CAMB code 23j based on 



the CMBFAST code 



24l | , where we include the modified Boltzmann equation of Eqs. ([T 



(IT9l)-(|2T]). Figure m shows the effects of m on the angular power spectrum with the following 
cosmological parameters: the baryon density parameter Vtji^ = 0.0223, the cold dark matter 
(CDM) density parameter Q^jjj^^h'^ = 0.104, K = 0, h = 0.73, and the reionization is 
neglected. We find that a decaying-A modifies the CMB power spectrum as follows: if m 
and/or is small, the amplitude of the power spectrum decreases. If we take larger values 
of m and/or i7^2) ^he first and third peaks of the power spectrum increase due to the large 
baryon density relative to the photon energy density. Furthermore, the CMB power spectrum 
shifts toward higher-/, because the photon last-scattering occurs at an earlier epoch as seen 
in Fig. [3] and Table [B We have found that the new term in the DACDM model increases 
the angular power spectrum at Z > 20 by about 10"^ percents. Therefore, the observational 
constraints are the same even if these terms are not included. 

The CMB angular power spectrum is rather sensitive to other cosmological parameters. 
For instance, baryon and CDM densities affect the amplitude of CMB anisotropy. Therefore, 
we need to carry out the MCMC approach 25| to constrain the possible parameters: i^^/i^, 
^1 reionization redshift z^g, the scalar spectral index n^, the amplitude of 
density fluctuation A^, and two parameters in DACDM (fi^2 ^^"^ note that we do 

not assume fiat universe. 

To start the MCMC calculations, we assume the priors on the cosmological parameters 
as follows: 0.5 <n^< 1.5, VL^^h'^ = 0.022 ±0.0022 at la C.L. (BBN prior), 0.01 < flcDAih^ < 
0.99, -0.3 < ^^/^o(= ~K/Hl) < 0.3, < fi^i < 1-0 , and 10 Gyr < to < 20 Gyr (age of the 
universe) . 
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FIG. 4: Comparison of the angular power spectrum in the decaying A model with the WMAP 
observation data. The solid line is the result of SACDM. The dashed, dot-dashed, and dotted lines 
are those of DACDM with {n^^^m) = (10-'^,0.5), (10-^0.8), and (10"^, 1.0), respectively. 



We constrain the relation between A2 (m) and other parameters from the recent CMB 



observations of WMAP 
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26j , BOOMERanG |27j , CBI |29j , and Acber 



28| . Our results 



are shown in Figs. [5] and [6l Figure [5] shows the constraints on the m — plane and our 
constraint is severer than that from the observed radiation temperature in FiglJl In our 
analysis, we obtain upper limits of m and f2^2 such as m < 4.2 and — ^-"^ ^ 10^^ at 
the 95.4 % confidence levels, respectively. Therefore, we cannot find the clear evidence of a 
decaying-A. Figure [6] shows the contours between fi^2 and other cosmological parameters 
{Vt^h? ,VL(^jjj^jh? ^rig, z^^Hq) at the 68.3% and 95.4% confidence levels. Parameters in the 
DACDM model have no degeneracy with other cosmological parameters; the parameters are 
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FIG. 5: Contours of the m — ^}j^2 plane from CMB. Soqlid lines indicate 68.3% and 95.4% confidence 
levels. Tlie dotted-line is the upper-hmit from Eqs. ([8]), ([9]), and ([TO]) with h = 0.73. 

independent on other parameters. Table HTl shows the comparison of cosmological parameters 
between the SACDM and DACDM model obtained from our MCMC analysis. As expected 
from Figs. [5] and [6] differences in the cosmological parameters are as small as some percents. 
However, the value of f2y\^ differs around 10 percents, which should be further constrained 
by future observations. 

IV. SUMMARY AND DISCUSSION 

We have investigated the possible difference in the thermal evolution of the universe 
with a decaying-A term as a function of the cosmic scale factor that reduces the photon 
energy density. Although the energy density of a A-term is increasing at the early era, the 
effects of the decaying A on the cosmic expansion rate can be ignored at the A dominate 
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FIG. 6: Constraints on logi7A2 against Qf,h'^ (left top panel), QcDAih'^ (right top panel), rig (left 
middle panel), Zre (right middle panel), fi^^^ (left bottom panel), and Hq (right bottom panel) 
from WMAP three year results. Two curves in each panel correspond to the 68.3% and 95.4% 
confidence levels, respectively. 



epoch. On the other hand, a decaying-A term has been insisted to alter the evolution of 



the photon temperature 
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14| . Depending on parameters in the DACDM model, the 



photon energy density could be lower or higher compared with that in the SACDM model at 
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2; > 0. However, the second case should not occur, because the photon temperature becomes 
negative at some epoch of z < 0. We estimate the epoch of the last-scattering surface by 
calculating the visibility function, and find that z^ec = 2040 for m = 1.2 and fi^g — 10~^ 
in the DACDM model, which indicate that the decoupling occurs earlier by Az^ec = 950 
compared to the case of the SACDM. 

We examine qualitatively the effects of a decaying-A term on CMB angular power spec- 
trum. We obtain the modified Boltzmann equation of photons in the DACDM model based 
on the cosmological perturbation theory and calculate the CMB angular power spectrum. 
We find that a decaying-A could alter the CMB angular power spectrum significantly due 
to the following reasons: large baryon energy density relative to the photon density causes 
to boost up the first and third peaks; the early photon decoupling shifts CMB spectrum to 
higher multipoles. 

Finally, using the Markov-Chain Monte Carlo analysis, we can put constraint on m, 
fi^2 5 ^^"i cosmological parameters. We obtain the upper limits of parameters in DACDM: 
m < 4.2 and Qj^2 < 1-7 x 10^^ 32] ■ In our analysis, the upper limit of m is close to 4. A 
decaying-A might affect primordial nucleosynthesis because the A-term evolves as radiation 
(photon, neutrino, and electron-positron) at the early universe. However, if m is large, A2 
tends to be small as shown in FigjSl As the result, A-term becomes smaller than radiation 
components. In fact, at T = 10^ K, the ratio of the A-term and radiation density is 10~^ 
for m = 4.2 and Qj^2 — 10^^*^ (upper limits from WMAP results). Therefore, we can say 
that effects of the early universe such as nucleosynthesis is small. Therefore, the effects of 
the A-term on the physical processes in the early universe, such as on the nuclosynthesis, is 
negligible 

Interestingly, there is no degeneracy between the two parameters in DACDM and other 
cosmological parameters. From these constraints, the contribution of a decaying-A term to 
the cosmic thermal evolution should be extremely small, since the best-fit values of m, 
are nearly zero. 

We assume a cosmological term as a function of scale factor for simplicity Even if we 
parameterize reasonably the evolution of the cosmological term or the equation of state of 
dark energy, our results would not change qualitatively. 

On the other hand, we find that the reionization occurs at z^e = 11 in the DACDM model, 
which suggests that a first object could be formed at around this epoch. We should note 
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TABLE II: Comparison of cosmological parameters between the SACDM and DACDM models 
obtained from Markov-Chain Monte Carlo analysis. 





DACDM 


SACDM 




< 1.7 X 10-^ 




m 


< 4.2 




logio (101%) 

Zre 
Us 

Ho 


n 0991 +0-0019 

U.U^Z1_Q QQ28 

ins+o-02i 

9 qoT+O.Oig 

10 6+'^-2 

945+0.050 
u.y*d_Q 031 

70.7 ± 19.7 


0.0223 ± 0.0007 
0.1037 ±0.0081 
3.156 ±0.056 

10.911? 
0.951 ±0.016 
71 ±3 


Qa 


n 7c;7+0.083 
•J- '^'-0.206 

249+0-261 


0.763 ± 0.034 
n oqq+o.oas 



that we assume that the reionization history can be described by step-function as discussed 
in Ref. [3]]. The next CMB satellite, Plank, is expected to determine detailed reionization 
history. Then a variable A-term model such as DACDM should be constrained further. 
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